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We study the spectrum of screening masses in the deconfined phase of (3+l)D SU(2) pure gauge theory near 
criticality and compare it with the spectrum of bound states in the broken symmetry phase of the 3D Ising model, 
which is related to the gauge theory by universality. 
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1. INTRODUCTION 



It is well known that both the 3D Ising model 



S 



E 



<n,m> 



and the 3D lattice regularized 



b A theory 



undergo a second order phase transition at /3 = j3 c 
from a phase where the Z(2) symmetry is broken 
((3 > (3 C ) to a phase when the Z(2) symmetry is 
restored (j3 < (3 C ). The order parameter of this 
transition is the expectation value of the funda- 
mental field, i.e. the magnetization. 

The spectrum of massive excitations in the bro- 
ken symmetry phase of 3D Ising model and </> 4 
theory can be determined from the analysis of 
wall-wall correlations of spin operators. By per- 
forming Monte Carlo simulations near /3 C , it has 
been possible to detect three + states M. Re- 
sults for their mass ratios agree in the two models: 
m2 /mi = 1-83(3) and m 3 /mi — 2.45(10) (mi is 
the fundamental mass). Since the 3D Ising model 
and the 3D (lattice) <fi 4 theory belong to the same 
universality class, this result shows that univer- 
sality is valid far beyond critical indices. 
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The first excitation (1712) is below the two par- 
ticle production threshold, it represents therefore 
a non-perturbative state, identified as a bound 
state of the fundamental excitation [g]. The sec- 
ond mass state (m.3) could be a fictitious one 
which mimics the two-particle cut. 

According to duality, the broken phase of the 
3D Ising model can be put in correspondence with 
the confined phase of the 3D Z(2) gauge theory. 
This implies that mass ratios in the glueball spec- 
trum in the confined phase of 3D Z(2) gauge the- 
ory should agree with mass ratios in the broken 
phase of 3D (spin) Ising model. This has been 
verified in Ref. using Monte Carlo results (|] 
for the + glueball masses in the 3D Z(2) gauge 
theory: m 2 /mi = 1.88(2) and m 3 /mi = 2.59(4). 

According to universality, the broken phase of 
3D Ising model corresponds to the deconfined 
phase of the (3+1 )D SU(2) pure gauge theory at 
finite temperature with the magnetization of 
the spin model mapped into the Polyakov loop, 
P(x,y,z) = TrH% =1 U4,(x,y,z,an4,). This im- 
plies that the same mass ratios observed in the 
broken phase of 3D Ising should appear in the 
spectrum of the slopes of the exponential decays 
of wall-wall correlations of Polyakov loop oper- 
ators, i.e. in the spectrum of screening masses, 
in the deconfined phase of (3+l)D SU(2) at fi- 
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nite temperature near criticality. In this work we 
present a numerical test of this prediction. 

2. NUMERICAL RESULTS 

A naive way to determine + screening masses 
is to study the connected wall-wall correlation 
function 

G(\ Z1 - z 2 \) = (P(zi)P(z 2 )} - (P( Zl )){P(z 2 )) 

Here P(z) = l/(N x N y )^ y P(x,y,z) is the 
wall-average of Polyakov loops, to is the low- 
est mass in the + channel, to* the first excited 
level in the + channel and the dots represent 
the contribution from higher mass excitations and 
from multi-particle cuts. Since (P(z)) ^ in 
the high temperature phase, in principle also a 
z-independent term in G(|.zi — z 2 \) of the form 
Co exp[— mN z ] should be included (see ||). How- 
ever, such term is evidently sub-dominant for 
to* < 2to, as we expect to be the case here. 

We performed Monte Carlo simulations using 
the overheat-bath updating algorithm j(| with 
Kennedy- Pendleton improvement 0. We stud- 
ied the naive wall- wall correlation at (3 — 2.33 
on a 18 2 x 36 x 4 lattice (statistics 1M) and at 
P = 2.36 on a 18 3 x 4 (statistics 1.5M). These 
values of (3 are slightly above the critical value 
(3 C = 2.29895(10) determined in Ref. §]. We de- 
termined the effective mass defined as 



G{z) =C(e 



) 



and found in both cases that data for m eff as a 
function of z reach a plateau value correspond- 
ing to the fundamental mass: m = 0.3660(20) at 
f3 = 2.33 and m = 0.5356(26) at (3 = 2.36. The 
deviations from the plateau at the smallest values 
of z can be attributed to lattice artifacts and to 
the possible effect of other physical states. In or- 
der to single out the latter effect, we rescaled the 
values of m otf and z by m = l/£ and £, respec- 
tively, and put together on the same plot data 
from (3 = 2.33 and (3 = 2.36. As shown in Fig. 0, 
data from the two different /3's fall on the same 
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Figure 1. Effective screening masses in physical 
units at j3 = 2.33 on a 18 2 x 36 x 4 lattice and 
13 = 2.36 on a 18 3 x 4 lattice. 

curve before reaching the plateau, indicating the 
physical origin of these effects. 

In order to detect this excited state, we adopted 
the so-called "variational" method ||,|l(|. The 
method includes the following steps: define a ba- 
sis of (wall-averaged) operators {O q }, compute 
the connected cross-correlation matrix among 
them 

C af3 (z) = (O a (z)Of,(0)) - {O a (z))(O p (fl)) , 
diagonalize C a ^(z) to obtain the eigenvalues 

Xi(z) cx e~ m ' z +e -™^ z -z) 

and finally extract the masses to^. In practical 
simulations not more than two (sometimes three) 
leading eigenvalues allow to extract a signal. The 
effectiveness of the method relies on the choice of 
a "good" set of operators; it is convenient to de- 
fine operators living on different length scales (for 
instance, by use of recursive definitions). More- 
over, by defining operators with non-trivial trans- 
formation under spatial rotation, it is possible to 
look for states with non-zero angular momentum. 

As a first set of operators we considered the one 
which can be built adapting to the present case 
the recursive procedure used in Ref. 0] : 



pW(x,y,z) = P(x,y,z) 
P {n+1) (x,y,z) = sign(w) 



(1 — w)\u\ + wy 
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u = j(p {n) (x-a,y,z) + P ( - n '>{x,y-a,z) 
+ P^(x + a,y,z) + P^(x,y + a,z^j , 

with w — 0.1 and y — (P(x,y, z)). We considered 
a set of six such (wall-averaged) operators corre- 
sponding to the smoothing steps 0, 3, 6, 9, 12 and 
15. The effective masses extracted by the varia- 
tional method applied to these set of operators 
are shown in Fig. |^ which shows a clear plateau 
at a value to(0 + ) = 0.3629(50), in agreement with 
the naive determination of the fundamental mass. 
Data coming from the next-to-leading eigenvalue 
and corresponding to the first excited level are 
quite noisy. Nevertheless, if we take as "plateau" 
value the first value of the effective mass which 
agrees at la level with the preceding and the 
following value, we get m*(0 + ) = 0.684(48) and 
a mass ratio m*(0 + )/m(0 + ) = 1.89(16) in good 
agreement with the prediction from universality. 
Further analysis has shown that this excited state 
is coupled with the operators with larger scale in 
the set we considered. Indeed, taking into ac- 
count only two operators (smoothing level and 
15) the same result is obtained. Instead, perform- 
ing the analysis taking only the first four opera- 
tors (smoothing levels 0, 3, 6, 9) gives m(0 + ) = 
0.3615(40) and m*(0+) = 0.848(28), with a mass 
ratio m*(0+)/m(0+) = 2.35(10) which is com- 
patible with the value of m 3 /mi in the 3D Ising 
model and could be the analog of the two-particle 
cut in the same model. 
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Figure 2. Lowest and first excited effective screen- 
ing masses in the + channel at (3 = 2.33 on a 
18 2 x 36 x 4 lattice from a set of six operators 
(smoothing levels 0, 3, 6, 9, 12 and 15). 



Finally, in order to look for states with non- 
zero momentum, we considered the following set 
of 2 + operators, inspired by Ref. H, 



x.y 



P(x, y, z)P{x + na, y, z) 



P(x,y,z)P(x,y + na,z) 



with n — 1,...,5. The variational method 
gave for the fundamental 2 + mass the value 
m(2+) = 1.23(14). The ratio of this value to 
the fundamental mass in the + channel (from 
the naive determination) m(0 + ) = 0.3660(20) 
gives m(2+)/m(0+) = 3.36(40). This result can 
be compared with the corresponding mass ra- 
tio in the glueball spectrum of the 3D gauge 
Ising model which, by duality, is connected to 
the spectrum of massive excitations in the bro- 
ken phase of the 3D (spin) Ising model. In the 
3D gauge model we have m(2 + )/m(0 + ) = 2.59(4) 
and m*(2+)/m(0+) = 3.23(7), so our result for 
the lowest screening mass in the 2 + channel seems 
to be in agreement with the first excited mass 
m*(2 + ), more than with the fundamental m(2 + ). 
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